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SECTION  I 
INTRODUCTION 


Considerable  thought  was  given  to  the  basic  approach  to  be  taken  in  the 
supersonic  analysis  of  interfering  axisynnnetric  bodies.  Initial  considera- 
tion was  given  to  a source  paneling  approach  since  this  theory  ideally  allows 
one  to  handle  nonaxisymmetric,  nonslender  bodies  of  any  configuration. 
Woodward's  theory  (Reference  1)  supposedly  is  valid  for  any  configuration; 
however,  only  source  panels  which  lie  parallel  to  the  free  stream  are  appli- 
cable to  his  theory.  Additionally,  it  appears  that  solutions  which  could  be 
developed  would  Involve  the  inversion  of  large  matrices.  Consequently,  this 
approach  was  abandoned  in  favor  of  a more  flexible  solution. 

The  second  approach  taken  in  the  analysis  was  the  utilization  of  linearly 
varying  line  source  distributions  placed  along  the  body  centerline  to  generate 
the  isolated  body  solution  and  supersonic  point  sources  strategically  placed 
along  the  image  line  as  determined  by  subsonic  theory  (References  2 and  3)  to 
counteract  the  presence  of  adjacent  bodies.  The  basic  theory  for  this  analy- 
sis is  presented  in  the  following  sections,  but  for  reasons  which  will  become 
apparent,  this  approach  fails  to  adequately  model  the  flow  field. 

The  third  approach  taken  in  the  analysis  is  similar  to  the  second  ap- 
proach except.  Instead  of  using  point  sources  to  counteract  the  interference 
effects,  linearly  varying  line  source  distributions  were  used.  The  results 
and  the  theory  for  this  technique  are  also  presented  in  the  following 
sections. 
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SECTION  II 


ISOLATED  BODY  SOLUTION 


The  concept  of  adding  solutions  to  generate  specified  body  shapes  in 
!•'  1 potential  flow  is  well  established  in  the  literature  even  for  linearized 

supersonic  flow  (References  4 and  5).  In  most  cases,  however,  only  the 
essential  features  of  the  derivation  are  presented.  For  purposes  of  this 
investigation,  certain  elements  of  the  derivation  are  expanded,  or  revised, 
for  the  interfering  body  cases.  In  order  to  provide  good  continuity  for  the 
whole  supersonic  analysis,  the  theoretical  derivation  for  an  isolated  thin 
supersonic  body  is  also  presented.  The  approach  as  presented  herein  for  a 
thin  isolated  body  closely  follows  that  of  Liepmann  and  Roshko  (Reference  4) . 

1.  Theory 

The  potential  for  tain  todies  in  supersonic  flow  is  governed  by  the  well' 
known  wave  equation 


3r2 


= 0 

9x2 


(1) 


It  has  been  shown  by  von  KirmAn  and  Moore  (Reference  5)  that  the  potential  ma’ 
be  written  as  an  integral  over  a distribution  of  sources  of  the  form 


<(>(x,r) 


X- 
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f(0  dg 

/(x-C)2  - A^r2 


(2) 


where  field  points  are  represented  by  (x,r)  and  5 is  the  coordinate  along  the 
source  distribution  as  shown  in  Figure  1.  Points  outside  the  zone  of  action 
or  outside  the  Mach  cone  produce  Imaginary  solutions  for  (Ji,  and  consequently 
only  those  points  which  lie  inside  the  Mach  cone  are  considered. 

It  is  convenient  to  make  a variable  transformation  of  the  form 


C = X - Ar  cosh  o , 


(3) 


Mach  cone  froir  orlpin 


Mach  cone  from  ^ 


Figure  1.  Line  Source  in  Supersonic  Flow 


in  which  case  the  potential  reverts  to 


In  Equation  (4),  f(?)  represents  the  function  which  governs  the  strength 
distribution  of  the  line  sources.  The  simplest  source  distribution  which 
produces  meaningful  results  is  (Reference  4) 


which  obviously  is  a linearly  varying  source  distribution  with  zero  strength 
at  the  nose.  For  a single  source,  this  solution  produces  a conical  body 
shape  with  the  flow  velocities  parallel  to  the  surface  everywhere  on  the  body 
Since  the  differential  equation  is  linear,  solutions  may  be  added  so  that  an 
arbitrarily  shaped  thin  axisymmetric  body  may  be  generated  by 


The  unknowns  in  Equation  (6)  are  the  a's  which  are  the  respective  slopes  of 
the  line  sources.  These  unknown  slopes  are  determined  by  matching  flow 
boundary  conditions  at  selected  points  along  the  body. 

Because  the.  flow  is  supersonic  and  sources  cannot  influence  points  up- 
stream of  their  own  Mach  cone,  a marching  technique  may  be  employed  in  which 
the  a’s  may  be  computed  in  a successive  fashion.  This  unique  feature  allows 
the  use  of  many  more  control  points  than  would  otherwise  be  practical  for 
matrix  inversion  techniques  using  the  same  amount  of  computer  time. 

Consider  the  case  shown  in  Figure  2. 


r 


Figure  2.  Schematic  of  the  Nose  of  a Body  in  Supersonic  Flow 


As  an  example,  consider  the  field  point  (x2*r2).  The  line  sources  which  in- 
fluence (x2,r2)  are  those  associated  with  a^^  and  a2.  The  line  source  of  slope 
a.^  begins  at  the  nose  or  at  and  the  line  source  of  slope  a2  begins  at  ^2* 
Note  from  Figure  2 that  the  locations  of  the  beginning  points  of  all  the  line 
sources  are  determined  from 


The  potential  at  any  field  point  affected  by  ai  and  a-2  is 


®1 


and  Equation  (4)  becomes 


/ 


aj^(x-Cj^  - Xr  cosh  a)  da 


o=cosh 


The  result  of  the  integration  is 


♦.j  ■ - [“»»■“  (-ir)-  ] 


The  potential  due  to  a2  is 


a2a-C2> 


2 J - x^r^ 


Using  the  variable  transformation.  Equation  (3),  the  resultant  potential  is 


[“»>’■' (tt)  - 


or,  in  general,  the  potential  due  to  any  jth  line  source  is 


.j  ■ - *j  ‘■‘■'j’  [ (^)  ■ '(j%)  ] 
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The  points  (x2»r2) , . . . (x^,rj^)  as  illustrated  in  Figure  2 are 

control  points  where  the  source  strengths  a]^,a2»...aj  are  determined  so 
that  the  velocity  vector  Is  parallel  to  the  local  body  slope  at  the  control 
points.  As  in  the  subsonic  case,  the  boundary  condition  is 


J^Vi 


(18) 


For  the  first  control  point  (xj^,rj^),  the  only  Influencing  line  source  is 
due  to  aj^.  Equation  (18)  then  becomes 


Rearranging  and  solving  for  a^^  yields 


For  the  second  point,  the  procedure  is  repeated  with  the  only  unknown  being 
B2  and  so  on  down  the  body.  In  general,  the  solution  for  the  ith  source 
(slope)  strength  is 

’•“i  4^) 

a ■ == (21) 


From  Equation  (21),  the  strength  of  all  the  N line  sources  are  computed  at 
N control  points. 

Once  all  the  strengths  are  known,  the  isolated  body  pressure  distribu- 
tion is  determined  using  Equations  (16)  and  (17)  as 


C = - 2u  + Xu^  - . 

P 


2.  Results  for  Body  Shape  A4V1P 

In  order  to  clearly  validate  the  isolated  body  solution  and  to  provide 
a check  case  for  the  computer  program,  a thin  axisymmetric  body  was  chosen 
for  which  pressure  data  was  available  in  the  low  supersonic  Mach  number  range. 
The  body  so  chosen  was  designated  as  the  A4V1P  (Reference  6).  Dimensions  for 
the  body  were  nondlmenslonalized  with  respect  to  the  M117  maximum  body  diame- 
ter in  order  to  provide  a standardized  means  of  comparison.  The  results  for 
the  Isolated  body  for  the  A4V1P  are  represented  in  Figures  3 and  4.  Note  that 
the  theoretical  results  are  well  within  experimental  error  when  compared  to 
wind  tunnel  results.  There  does  seem  to  be  a slight  tinderprediction  by  the 
theory  in  the  region  of  the  tail  which  is  probably  attributable  to  buildup  of 
the  boundary  layer  in  this  region.  As  a general  rule,  however,  there  is  ex- 
cellent agreement  between  theory  and  experiment  for  the  Isolated  A4V1P  body. 


□ Experimental 


Figure  4 


Comparison  of  Theoretical  and  Experimental  Pressure 
Distribution  for  the  A4V1P  Missile  at  a=0.0  and  M=1.56 


3.  Results  for  the  M117  Body  Shape 


The  second  body  chosen  for  the  analysis  was  the  M117  body  shape.  This 
body  was  chosen  in  order  to  illustrate  the  effectiveness  of  the  theory  for 
bodies  with  blunter  noses,  smaller  L/D  ratio  bodies,  and  bodies  with  discon- 
tinuities in  the  surface  slope.  As  for  the  previous  body,  linear  dimensions 
were  nondimensionalized  with  respect  to  the  M117  maximum  body  diameter  of 
16  Inches. 

For  the  M117,  the  Mach  number  range  was  limited  to  a maximum  of  1.5.  The 
results  are  presented  in  Figures  6 through  8 . Obviously,  the  M117  is  not  a 
slender  body,  and  at  these  Mach  numbers  the  nose  shock  is  detached  from  the 
body.  However,  there  is  nothing  in  the  Isolated  theory  to  account  for  the 
presence  of  shock  waves.  Figure  6 is  a comparison  of  the  theoretical  and 
e:q)erlmental  results  for  the  M117  at  a Mach  number  of  1.1  (References  7 and 
8).  Since  this  clearly  is  transonic,  mixed  subsonic  and  supersonic  flows, 
theoretical  results  might  be  questionable;  however,  the  theory  agrees  remark- 
ably well  with  wind  tunnel  data  (Reference  7).  No  pressure  data  was  available 
for  the  M117  at  Mach  numbers  greater  than  1.1. 

Even  though  shock  waves  are  not  modeled  in  the  analysis,  there  is  a limit 
on  the  Mach  number  at  which  the  M117  can  be  theoretically  analyzed  with  line 
sources  placed  along  the  body  centerline.  The  slope  or  nose  angle  of  the  M117 
is  approximately  41.41  degrees  measured  from  the  (horizontal)  centerline  axis. 
Consequently,  any  Mach  number  with  associated  Mach  angle  B which  is  less  than 
41.41  degrees  causes  the  wave  angle  to  intersect  the  body  as  illustrated  in 
Figure  5,  below.  The  Mach  number  associated  with  a wave  angle  of  41.41“  is 


Figure  S.  Nose  and  Centerbody  Section  of  the  M117  Store 
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1.51.  For  Mach  numbers  greater  than  1.51,  the  associated  Mach  angle  is  less 
than  the  nose  angle;  consequently,  the  theoretical  analysis  becomes  inad- 
equate. For  this  reason,  only  Mach  numbers  up  to  1.5  were  considered  for 
the  M117  store. 


In  summary,  the  technique  of  using  line  source  distributions  to  model  a 
thin  body  in  supersonic  flow  is  well  established.  Additionally,  the  theory  as 
herein  applied  adequately  predicts  the  basic  axial  pressure  distribution  for 
thin  bodies  and  appears  adequate  for  bodies  such  as  the  M117.  However,  the 
applicable  Mach  number  range  is  limited  by  the  nose  shape,  and  for  nose  shapes 
blunter  than  the  M117,  the  Mach  number  range  is  severely  restricted. 


SECTION  III 


TWO-BODY  SOLUTION 


There  were  two  attempts  to  model  the  two-body  interference  effects  in 
supersonic  flow,  the  first  using  a point  source  image  solution  and  the 
second  using  a linearly  distributed  line  source  solution.  Although  the  point 
source  solution  did  not  adequately  model  the  flow  field,  the  basic  theory  and 
associated  results  are  presented,  and  the  reasons  for  its  failure  are  also 
discussed. 

1.  Point  Source  Image  Solution 

For  the  beginning  of  the  interference  problem,  it  is  assumed  that  the 
isolated  body  solution  using  linearly  distributed  line  sources  is  completed. 
For  the  interference  effects,  point  sources  are  placed  along  an  image  line  in 
each  body  to  counteract  the  crossflow  velocities  and  maintain  the  tangency 
boundary  conditions  at  a set  of  control  points  on  the  surface  of  each  body 
(see  References  2 and  3).  The  two  bodies  are  arranged  as  shown  in  Figure  9. 
The  image  line  is  defined  as  in  subsonic  theory  as 


6 


(1.0  - /1.0-(2r/C^^)2 


(23) 


In  order  to  preserve  the  short  computational  time,  a scheme  was  devised 
so  that  the  marching  technique  could  still  be  employed  for  the  image  system 
solution  as  was  done  for  the  isolated  body  solution.  Referring  to-Figure  9, 
the  computational  steps  proceeded  in  the  following  manner. 

Geometry  - (1)  The  point  labeled  "0"  on  top  of  body  1 (the  lower  body) 
was  determined.  This  is  the  point  on  body  1 where  the  influence  of  body  2 
is  first  felt.  This  point  is  the  Intersection  of  the  downstream  running 
Mach  wave  from  the  nose  of  body  2 and  the  top  of  body  1.  (2)  Next,  the  point 

labeled  "1"  on  the  image  line  of  body  1 is  determined  as  the  intersection  of 
the  upstream  running  Mach  line  from  the  zero  point  (on  top  of  body  1 just 
determined)  and  the  image  line,  as  defined  in  Equation  (23),  above.  This  is 
the  point  at  which  the  first  point  source  is  to  be  placed.  (3)  The  next 
step  is  to  determine  point  1 on  the  centerline  of  body  1 and  point  zero  on 
the  bottom  of  body  1.  The  point  1 on  the  centerline  of  body  1 is  the  loca- 
tion of  the  first  sink  to  approximately  counteract  the  effects  of  the  first 
source  on  the  image  line.  (4)  Point  1 on  the  top  of  body  1 is  located  at 
the  same  axial  location  as  point  zero  on  the  bottom  of  body  1.  (5)  With 

point  1 on  top  of  body  1 determined,  point  2 on  the  image  line  is  then  de- 
termined and  the  whole  process  is  repeated.  As  each  point  in  body  1 is  de- 
termined, the  corresponding  point  in  body  2 is  also  stored  in  computer 
memory. 


Figure  9.  Schematic  of  Two-Body  Interference  Geometry  Using  Point  Sources 


Computational  Procedure  - (1)  The  first  control  point  on  top  of  body  1 
is  selected.  The  induced  velocities  are  determined  from  the  line  sources  in 
body  2,  any  point  sources  in  body  2,  the  line  sources  in  body  1,  and  any 
point  sources  in  body  1.  From  geometry  (Figure  9),  it  is  observed  that  the 
only  unknown  point  source  which  can  influence  the  first  control  point  on  top 
of  body  1 is  the  point  source  located  on  the  image  line  of  body  1;  the 
strength  of  this  point  source  is  determined  to  maintain  flow  tangency  at  the 
top  control  point.  (2)  The  first  control  point  on  the  bottom  of  body  1 is 
next  selected.  The  same  procedure  of  determining  induced  velocities  from  all 
upstream  line  and  point  sources  is  followed  for  this  control  point  on  the 
bottom  of  body  1.  The  only  unknown  point  source  strength  in  the  upstream 
running  Mach  cone  from  the  point  1 (bottom)  control  point  is  the  one  located 
at  the  point  1 centerline  location.  Consequently,  the  strength  of  this  source 
is  deterniined  to  maintain  flow  tangency  on  the  bottom  of  body  1 at  the  first 
control  point.  (3)  Next,  the  second  control  point  on  top  of  body  1 is 
selected,  and  the  strength  at  the  point  2 image  location  is  determined.  In 
this  manner  alternately  selecting  control  points  on  top  and  then  on  bottom  of 
body  1,  the  strengths  of  the  point  sources  located  along  the  image  line  and 
along  the  centerline  are  respectively  determined  which  forces  the  flow  to  be 
tangent  to  the  body  both  on  top  and  on  bottom  of  each  body  at  each  control 
point . 

Results  of  Point  Source  Solution  - The  resulting  solution  from  the  point 
source  analogy  hinges  on  the  basic  formulation  of  the  induced  velocity  from 
a supersonic  point  source.  The  potential  for  a supersonic  point  source  is 

4>  = - m/h  (24) 


where  m is  the  source  strength  and  h is  the  hyperbolic  radius  defined  by 


h = /(x-x  )2  -X2[(y-y  )2  + (z-z  )2]  . (25) 

The  resulting  induced  velocities  are 


= = m(x-x^)/h3  , 

(26) 

= 11=  _ mX2(y-y^)/h3  , 

(27) 

= = - mA^(z-z  )/h3  . 

oZ  O 

(28) 

In  the  above  equation,  x , y , z is  the  location  of  the  source  and  x,  y,  z 
is  the  field  point.  o o o 
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Two  basic  observations  should  be  noted  from  Equations  (26)  through  (28): 
(a)  the  induced  velocity  is  dependent  on  how  close  the  field  point  is  to  the 
source  location  and  (b)  the  induced  velocity  is  strongly  dependent  on  how 
close  the  field  point  is  to  the  Mach  line.  As  a matter  of  fact,  as  the  field 
point  gets  close  to  the  Mach  line,  h approaches  zero,  and  the  induced  velocity 
approaches  infinity.  As  clearly  illustrated  in  Figure  9,  upstream  running 
Mach  lines  from  the  control  points  come  very  close  to  sources  in  body  2.  Con- 
sequently, large  induced  velocities  at  some  of  the  control  points  were  en- 
countered with  a subsequent  large  source  strength  to  offset  the  induced 
velocities.  Whenever  one  exorbitant  source  strength  was  determined,  each 
successive  strength  became  larger,  thereby  causing  the  solution  to  blow  up. 
Increasing  the  number  of  images  and  control  points  only  compounds  the  problem, 
and  reducing  the  number  of  images  causes  exorbitant  strengths  due  to  the  very 
wide  spacing  of  the  image  system.  Consequently,  this  approach  was  abandoned. 


2.  Line  Source  Image  Solution 


I As  in  the  point  source  approach,  it  is  assumed  that  the  isolated  body 

; solution  has  already  been  completed  by  placing  linearly  varying  sources  along 

the  centerline  of  each  body.  The  point  on  body  1 where  the  nose  of  body  2 is 
first  felt  is  determined  as  before  by  finding  the  intersection  of  the  Mach 
line  from  the  nose  of  body  2 and  the  top  of  body  1.  The  first  control  point 
for  the  interference  system  is  then  selected  as  the  next  evenly  spaced  axial 
, location.  For  example,  in  Figure  10,  the  standard  evenly  spaced  x locations 

j are  0.2,  0.4,  0.6,  0.8,  and  so  on  down  the  body.  The  axial  location  where 

the  Interference  from  the  nose  of  body  2 is  first  felt  on  body  1 is  approx- 
j imately  0.54  (Figure  10),  designated  point  "0"  on  top  of  body  1.  The  next 

I evenly  spaced  axial  location  is  0.6,  which  is  designated  as  point  1 on  the 

I top  of  body  1 and  point  0 on  the  bottom.  Evenly  spaced  axial  locations,  0.8, 

1.0,  1.2,  ...,  are  then  selected  to  define  other  downstream  control  point 
I locations  on  each  body. 
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Note  from  Figure  10  that  the  location  of  the  first  image  line  source  is 
the  intersection  of  the  upstream  running  Mach  line  from  point  0 on  top  of 
body  1 and  the  upstream  running  Mach  line  from  point  0 on  the  bottom  of  body 
1.  In  order  to  more  clearly  illustrate  the  geometry  associated  with  the  com- 
putational system.  Figure  11  is  a schematic  of  the  nose  section  of  one  of  the 
interfering  bodies  as  illustrated  on  Figure  10.  The  axial  distance  to  the 
beginning  point  of  the  image  line  is  designated  as  Cm  and  from  Figure  11, 


t^k-1  - - \^^top 


(29a) 


where  is  the  radial  distance  from  the  body  centerline  to  the  kth  image 
point.  Note  from  Figure  11  that 


^^-l^bot  ^^^top 

and 


^’^k-l^bot  ^*^k^top’ 


From  Equations  (29)  and  (30),  then 


(30a) 

(30b) 


(31) 


Each  line  source  extends  from  its  point  of  origin,  as  defined  by  Equations 
(29)  and  (31) , to  infinity  and  is  parallel  to  the  x axis  as  shown  in  Figure 
10.  As  the  origin  of  each  line  source  location  is  determined  for  body  1, 
corresponding  locations  in  body  2 are  also  stored  in  computer  memory.  Note 
that  the  location  of  the  image  line  is  determined  by  the  Ax  interval  chosen 
for  the  computations  and  not  from  a predetermined  equation.  Actually,  this 
poses  no  problem  because  velocities  induced  from  linearly  varying  line 
sources  are  constant  along  any  Mach  line  emanating  from  the  line  source,  i.e., 
(x/Ar)  = constant. 

The  locations  of  the  beginning  point  of  the  centerline  line  sinks 
are  determined  from  Figure  11  as 


^^c^^k  " *k  “ ^ ’^k  • 

These  sinks  along  the  centerline  are  used  to  match  boundary  conditions  on 
the  bottom  of  body  1 due  to  the  presence  of  body  2. 

The  equations  which  govern  the  induced  velocities  from  these  line 
sources  are  basically  Equations  (16)  and  (17)  which  are  slightly  modified 
as 


and 


cosh~^ 


(33) 


V 


i 


(34) 
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where  rp  is  the  radial  distance  between  the  control  point  and  the  begin- 
ning of  the  line  source  (i  = control  point) . It  should  be  noted  from  the 
geometry,  as  illustrated  on  Figure  10,  that  a marching  technique  can  be 
utilized  to  find  the  unknown  strengths  of  the  image  sources  to  match  boundary 
conditions  on  top  of  body  1 and  the  unknown  sink  strengths  placed  along  the 
body  centerline  to  match  boundary  conditions  on  the  bottom  of  body  1. 

One  additional  restraint  is  placed  on  the  marching  scheme  from  Figure  10; 
consider,  as  an  example,  control  point  designated  number  3 on  top  of  body  1. 
Note  that  the  upstream  running  Mach  line  from  this  control  point  encounters 
the  image  line  source  number  3 in  both  body  1 and  body  2.  Since  this  double 
encounter  may  or  may  not  occur  for  any  one  particular  control  point,  pro- 
visions must  be  made  in  the  computational  procedure  to  include  the  double 
Induced  velocity  for  the  affected  control  points.  The  technique  used  in  this 
analysis  is  outlined  below.  (1)  An  array  is  set  up  in  the  computer  which 
contains  the  strength  of  the  image  line  sources.  This  array  is  initially  set 
up  with  zeros  stored  in  all  of  its  locations.  (2)  As  each  line  strength  is 
determined,  the  corresponding  zero  is  replaced  by  the  computed  line  source 
strength.  (3)  Control  points  as  previously  determined  are  successively  chosen 
in  a downstream  marching  fashion.  Consider  as  an  example  the  kth  control 
point.  The  Induced  velocity  from  all  line  sources  in  the  upstream  running 
Mach  cone  is  determined  with  the  only  unknown  (appearing  in  both  bodies)  being 


the  kth  line  source  strength  a. 


m^- 


The  boundary  condition  then  becomes 


f 

I y, 


V,  + V 
k r 

1 + u,  + V 
k X 


(slope) 


(35) 


where  and  are  the  accumulated  sum  of  the  radial  and  axial  velocity  com- 
ponents of  all  known  line  sources  in  the  upstream  Mach  cone  of  the  kth  con- 
trol point. 


From  Equations  (33)  and  (34)  and  taking  into  account  the  double  encounter 
effect,  the  axial  and  radial  Induced  velocities  from  the  kth  line  source 
become,  respectively. 
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From  Equations  (35),  (36)  and  (37),  the  line  strength  is  isolated  as 
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S,  (1.0  + V )-V 
k X r 

"he  Xk,  + S,  k., 
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In  the  above  equation. 


/G^S-1  - - 1 

cosh“^  G^  + cosh“^ 

[(x-C)/(Xr)]Body  1 
[(x-C)/(Xr)lBody  2 


where 


and 


k.  = 


k„  = 


G,  = 


(39] 

(4o; 


(A2) 


For  control  points  where  the  double  encounter  is  not  a problem,  it  may  be 
shown  that 


S,  (1.0+V  ) - V 
k X r 

"he  X/Gj^2-1  + cosh“l  Gj^ 


(43) 


Centerline  strengths  are  determined  in  a similar  manner  to  match  boundary 
conditions  on  the  bottom  of  body  1.  The  results  of  the  computations  for 
the  interference  solution  are  presented  in  the  next  section. 


3.  Results  for  Body  Shape  A4V1P 

For  the  A4V1P,  it  was  necessary  to  assume  some  spatial  centerline 
distance  for  the  two  interfering  bodies . Since  there  was  no  experimental 
data  for  comparison  purposes,  the  nondiroensional  centerline  distance  was 
chosen  to  be  1.0.  This  choice  provided  good  Interference  effects  even  on 
the  nose  of  the  two  bodies. 


The  results  for  the  interference  pressure  distribution  is  shown  on 
Figures  12  through  16.  The  angular  locations  refer  to  body  1 (bottom 
body)  with  boundary  conditions  being  matched  exactly  on  top  of  body  1 
(90°)  and  on  bottom  of  body  1 (-90°). 
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To  lend  credence  to  the  formulation  of  the  problem  and  the  resulting 
solution,  the  -90°  case,  bottom  of  body  1 on  Figure  14  or  15,  is  very  nearly 
the  same  as  the  free  stream  data  as  Illustrated  in  Figure  3.  As  the  two 
bodies  interfere  with  each  other,  the  pressure  on  top  of  body  1 begins  to 
increase,  meaning  that  the  flow  is  slowing  do\«m  or  approaching  M=1.0.  One 
would  expect  this  result  since  the  flow  is  converging  between  the  bodies. 

Note  also  that  the  resulting  pressure  distribution  on  the  body  produces  a 
very  large  negative  pitching  moment  which  is  in  agreement  with  experimental 
results  for  stores  of  similar  shape.  However,  since  no  actual  interference 
data  for  the  A4V1P  is  available,  comparison  of  the  magnitude  of  the  numbers 
for  the  Interference  case  is  not  possible. 

It  has  previously  been  stated  that  the  radial  position  of  the  image  line 
has  no  influence  on  the  induced  velocities.  Also,  from  the  derivation  of  the 
interference  equations,  the  radial  position  of  the  image  line  is  determined 
by  the  Ax  interval  chosen  for  the  analysis.  In  order  to  clearly  illustrate 
this  point.  Figure  14  is  the  interference  results  using  60  line  sources  and 
Figure  15  is  the  interference  results  using  120  line  sources.  The  difference 
in  the  pressure  distributions  from  these  two  runs  are  of  minor  magnitudes 
as  expected. 

Since  the  A4V1P  is  indeed  a slender  body,  the  theoretical  results  were 
well  behaved  with  pressure  distributions  much  as  expected.  Because  the  M117 
is  not  a slender  body  and  because  of  the  maximum  Mach  number  limit  encountered 
on  the  nose  for  the  Isolated  body  solution  as  previously  discussed,  one  might 
expect  additional  problems  for  interfering  M117  stores. 

4.  Results  for  the  M117  Body  Shape 

The  problem  of  the  limited  Mach  number  range  for  the  M117  associated  with 
the  nose  geometry  of  course  still  existed  for  the  interference  cases.  In 
addition  to  the  nose  angle  problem,  the  computations  for  the  flow  over  the 
sharp  shoulder  proved  to  be  somewhat  of  a problem  for  the  interfering  case. 

In  the  real  situation,  an  expansion  fan  must  exist  at  the  shoulder  if  the 
flow  is  to  remain  supersonic.  However,  the  free  stream  flow  must  possess  a 
sufficiently  high  Mach  number  if  it  is  to  make  the  turn  without  separating  or 
shocking  down.  For  the  isolated  body  (M117),  this  problem  did  not  emerge 
from  the  solution  as  is  shown  in  Figures  5 through  7;  however,  for  the 
two-body  case,  the  mathematical  analysis  does  Indeed  become  unstable  when 
computations  are  attempted  downstream  of  the  shoulder. 

In  Figures  17  through  19,  the  results  for  the  interference  calculations 
for  two  adjacent  M117  stores  are  presented.  The  centerline  distance  is  1.042 
diameters  as  is  approximately  the  case  when  the  stores  are  mounted  on  the 
triple  ejector  rack  (TER) . Note  that  for  these  plots  the  scale  on  the  pres- 
sure distribution  is  an  order  of  magnitude  smaller  than  on  the  A4V1P  re- 
sults. Note  also  that  the  computations  for  Cp  were  limited  to  ±2.0  since 
pressure  coefficients  outside  this  range  are  completely  unrealistic. 

At  very  low  supersonic  Mach  numbers  (Figure  17),  the  computational 
technique  becomes  unstable  aft  of  the  sharp  shoulder.  From  physical 


principles  of  the  flow  between  adjacent  bodies,  some  pressure  recovery 
mechanism  such  as  a standing  normal  shock  is  necessary  for  low  supersonic 
Mach  numbers.  Since  this  mechanism  is  not  incorporated  in  this  analysis, 
the  solution  becomes  unstable.  This  instability  is  also  seen  in  Figures 
18  and  19.  Note,  however,  as  the  free  stream  Mach  number  increases  (M=1.5), 
the  point  on  the  body  where  the  solution  becomes  unstable  moves  downstream 
of  the  shoulder. 

Specific  conclusions  regarding  the  agreement  of  the  theory  with  experi- 
mental dataware  not  possible  since  no  data  for  a two-store  M117  configuration 
was  available. 

5.  Surface  Velocities  on  the  M117 

Another  Interesting  aspect  of  these  theoretical  results  for  the  M117 
emerged  from  the  computations.  In  the  analysis  it  is  assumed  that  small 
perturbation  theory  is  applicable.  Therefore,  the  velocity  Induced  from  the 
line  sources  must  be  much  less  than  the  free  stream  velocity.  From  Equation 
(16)  it  is  noted  that,  for  a linearly  varying  line  source,  the  axial  component 
of  the  velocity  is  negative.  Consequently,  at  small  Mach  numbers,  i.e., 

M=l.l,  even  small  perturbations  may  generate  subsonic  surface  velocities.  If 
the  perturbation  velocities  are  to  maintain  supersonic  flow  on  the  body  sur- 
face, then 


(V  + u')2  + (v')2  > a2  (44) 

00  oo 

Nondimensionalizing  with  respect  to  yields 

(1  + u)2  + v2  > 1/m2  <^5) 

00 

The  left-hand  side  of  Equation  (45)  is  the  nondimens ional  surface  velocity, 
or 


V = /(I  + u)^  + . 

s 


(46) 


For  supersonic  flow  on  the  body  surface,  then 


As  an  example,  for  M=l,l,  Equation  (47)  becomes 

V > .9091  . 


(47) 


(48) 
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Figure  20  Is  a plot  of  the  nondimenslonallzed  surface  velocity,  Vg,  for  the 
M117  at  Mach  numbers  of  1.1  and  1.3.  This  figure  and  Equation  (47)  indicate 
that  subsonic  flow  must  exist  over  certain  portions  of  the  nose  if  the  body 
boundary  conditions  are  to  be  met.  This  simple  analysis  also  theoretically 
predicts  the  sonic  point  on  the  nose  and  the  subsequent  region  of  supersonic 
flow. 


The  predicted  areas  of  subsonic  and  supersonic  flow  appear  realistic; 
however,  It  should  be  stressed  here  that  the  authors  do  not  Imply  that  the 
results  are  rigorous  and  certainly  many  other  factors  should  be  included  in 
the  analysis.  If,  however,  the  results  were  taken  literally,  then  the  nose 
section  of  the  M117  should  be  treated  with  a subsonic  analysis  and  the  aft 
section  with  supersonic  analysis.  Future  work  in  this  area  could  result  in 
one  program  that  could  perhaps  handle  both  situations  or,  in  essence,  pro- 
vide a simplified  transonic  solution. 
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Figure  20.  Surface  Velocity  Distributions  on  the  Nose  Section 
of  the  M117  at  Macli  Nur.bers  of  1.1  and  1.3 
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SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  technique  of  using  a system  of  linearly  varying  line  sources  to 
generate  a specified  axisymmetric  body  shape  ^n  supersonic  flow  has  proven 
very  successful  (References  4 and  5).  Pressure  distributions  obtained  from 
this  theory  are  quite  good  and  are  within  experimental  error  for  thin  bodies. 

Computational  problems  arise  in  the  nose  region  of  thicker  bodies  such 
as  the  M117  store.  If  the  initial  slope  of  the  nose  is  larger  than  the  Mach 
angle,  centerline  source  filaments  cannot  be  used  to  model  the  body  shape 
because  portions  of  the  nose  lie  outside  of  the  zone  of  influence  of  the 
filaments.  This  is  indeed  the  case  for  the  M117  at  Mach  numbers  greater  than 
about  1.5.  Consequently,  the  M117  analysis  was  limited  to  Mach  numbers  below 
1.5. 


As  a basic  ground  rule  in  developing  an  interference  theory  for  super- 
sonic flow,  it  was  decided  that  the  downstream  marching  idea  should  be  main- 
tained since  other  schemes  using  matrix  Inversion  would  require  longer  com- 
putational times  and  increasing  storage.  Two  approaches  for  the  supersonic 
interference  problems  were  considered.  The  first  approach  used  supersonic 
point  sources,  placed  along  an  image  line  and  along  the  body  centerline,  to 
account  for  interference  effects.  Severe  computational  problems  were  en- 
countered with  this  approach  due  to  a mathematical  singularity  in  the  Induced 
velocity  at  control  points  or  field  points  lying  on  or  near  a Mach  line.  It 
is  difficult  to  generate  two  Interfering  bodies  in  which  a control  point  does 
not  lie  near  a Mach  line.  Consequently,  it  is  concluded  that  this  approach 
is  inadequate  in  modeling  the  flow  field  in  supersonic  flow  for  two  inter- 
fering bodies. 

The  second  approach  utilized  linearly  varying  supersonic  line  sources 
for  the  image  system.  The  unique  features  of  this  approach  are:  (1)  no 

singularities  exist  in  the  flow  field  even  on  the  Mach  cone,  and  (2)  induced 
velocities  are  constant  along  Mach  lines  emanating  from  the  line  source 
itself.  The  marching  technique  is  maintained,  and  no  assumptions  concerning 
the  location  of  the  image  line  are  necessary.  The  results  of  these  computa- 
tions for  thin  bodies  appear  realistic,  but  no  data  was  available  for  com- 
parison. For  bodies  such  as  the  M117,  there  are  computational  problems  when 
low  Mach  number  flow  expanded  around  sharp  corners  or  downstream  of  sharp 
comers . 

From  the  overall  results  thus  far  obtained,  several  interesting  concepts 
or  recommendations  have  emerged: 

1.  The  idea  of  using  supersonic  line  sources  and  perhaps  line  doublets 
to  model  interfering  bodies  appears  promising. 
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2.  Experimental  pressure  data  for  two  interfering  bodies  should  be 
obtained  to  verify  the  analysis. 

3.  Extension  of  the  theory  to  include  some  boundary  layer  effects  is 
necessary. 

4.  Theoretical  approaches  for  determining  the  Mach  number  distribution 
over  the  nose  section  of  blunt  interfering  bodies  must  be  developed. 

5.  The  potential  theory  technique  for  identifying  subsonic  and  super- 
sonic flow  regions  should  be  pursued  and  more  rigorously  founded  in  theory. 

6.  A complete  program  for  handling  simultaneously  both  subsonic  and 
supersonic  flow  regimes  on  a single  body  needs  to  be  developed  as  an  initial 
effort  in  the  solutions  of  the  transonic  flow  field. 
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i LIST  OF  VARIABLES 
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Variable 

Name 

A 

ACLl 

ACL2 

AML 

AM2 


Definition 

Slope  of  linearly  varying  line  sources  for  the  Isolated  body. 

Slope  of  linearly  varying  line  sources  for  the  centerline 
image  strengths  in  body  1. 

Slope  of  linearly  varying  line  sources  for  the  centerline 
image  strengths  in  body  2. 

Slope  of  linearly  varying  line  sources  for  the  image  strengths 
in  body  1. 

Slope  of  linearly  varying  line  sources  for  the  image  strengths 
in  body  2. 


CLD 


Nondlmenslonal  distance  between  body  centerlines  for  two  inter- 
fering bodies. 


Pressure  coefficient. 


INPUT  variable.  If  IPLOT  is  input  as  0,  no  plots  are  obtained. 
If  it  is  input  as  1,  plots  are  obtained. 

Axial  coordinate  nondlmenslonal  distance  from  the  nose  of  the 
body  to  the  first  shoulder  (end  of  ogive  nose  section) . 

Axial  nondlmenslonal  distance  from  nose  of  body  to  second 
shoulder  (end  of  cylindrical  section) . 

Axial  nondlmenslonal  coordinate  length  of  body. 


Mach  number. 


Number  of  centerline  sources  to  be  placed  in  the  body. 


Number  of  axial  points  on  the  body  at  which  the  pressure  is 
to  be  computed. 

Number  of  azimuthal  locations  around  the  body  at  which  pressure 
confutations  are  to  be  made.  The  angular  increments  are 
every  45 

Local  radius  of  body. 
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LIST  OF  VARIABLES  (CONCLUDED) 


Variable 

Name 


Definition 


Radial  distance  from  the  coordinate  x-axis  to  the  centerline 
images  in  body  1 (RCLl  = 0.0). 

Radial  distance  from  coordinate  x-axis  to  centerline  image 
sources  of  body  2 (RCL2  = CLD) 

Radial  distance  from  x-axis  to  image  line  of  body  1. 

Radial  distance  from  the  coordinate  x-axis  to  the  image 
line  in  body  2. 

Maximum  body  radius  in  nondimensional  units . 

Local  slope  of  body. 

Axial  distance  from  coordinate  axis  to  beginning  point  of 
centerline  line  sources. 

Tail  angle  of  M117  in  degrees. 

Tangent  of  tall  angle  on  the  M117. 

Azimuthal  angle  around  body  at  which  the  first  set  of  pressures 
are  to  be  computed. 


Local  radial  nondimensional  velocity  component. 

Local  axial  nondimensional  velocity  component. 

Axial  distance  along  the  body. 

Axial  distance  to  beginning  point  of  centerline  Images  in 
body  1 or  body  2. 

Axial  distance  from  coordinate  axis  to  beginning  point  of  image 
line  sources. 


INTRODUCTION 


A copy  of  the  computer  program  used  In  the  interference  analysis  is 
included  in  this  appendix.  The  case  as  illustrated  in  the  program  listing 
is  for  the  M117  store  at  Mach  1.3.  A brief  discussion  of  each  subroutine 
is  included  below  with  required  inputs.  A listing  of  the  major  variables 
in  the  program  is  also  included. 

MAIN 


The  executive  part  of  the  program  calls  the  required  subroutines  which 
perform  various  tasks  to  determine,  ultimately,  the  pressure  distribution 
on  two  Interfering  bodies.  The  general  sequence  of  events  which  MAIN  per- 
forms is  as  follows:  (1)  Sets  up  the  geometry  for  a single  body  solution. 

(2)  Generates  the  single  body  line  source  solution.  (3)  Sets  up  locations 
of  the  image  sources,  along  image  line  and  body  centerline  for  the  two-body 
interfering  case.  (4)  Determines  the  coordinates  of  the  point  where  body  2 
first  interferes  with  body  1.  (5)  Generates  the  interfering  body  solution 

for  the  image  source  strengths  and  the  image  centerline  source  strengths. 

(6)  Determines  the  pressure  distribution  on  the  body  at  selected  coordinate^ 
locations.  (7)  Plots  desired  results. 

The  inputs  to  MAIN  are  in  the  form  of  data  cards  as  shown  on  the  next 
page.  This  information  is  passed  to  the  various  subroutines  through  labeled 
COMMON  blocks  as  listed  at  the  beginning  of  MAIN. 

SUBROUTINE  GEOMB  (RAD,  X,  SLOPE,  NS) 

The  GEOMT-  subroutine  determines  the  local  body  radius  and  body  slope 
at  discrete  axial  locations.  The  axial  locations  enter  the  subroutine 
through  the  CALL  argument  list  as  the  X array.  The  array  which  contains 
corresponding  body  radii  is  the  RAD  array,  and  the  SLOPE  array  contains  the 
corresponding  body  slopes.  These  arrays  are  transferred  back  through  the 
CALL  argument  list.  Other  necessary  input  parameters  for  GEOMB  enter  through 
the  COMMON  statement.  The  variable  NS  in  the  CALL  argument  list  is  the 
number  of  axial  locations  contained  in  the  X array.  Other  variables  have 
already  been  identified  in  the  main  program. 

SUBROUTINE  RHP  (K,  X,  R,  XI,  RI,  THA,  ACC) 

The  RHP  subroutine  determines  how  many  of  the  line  sources  in  body  1 or 
2 Influence  any  particular  control  point  under  consideration.  This  number 
is  transferred  back  through  the  CALL  argument  list  as  K.  Other  arrays  which 
are  transferred  to  RHP  through  the  CALL  are  X,  R,  XI,  and  RI.  These  arrays 
are  generated  elsewhere  in  the  program.  The  variable  THA  in  the  CALL  argu- 
ment is  the  azimuthal  angle  around  the  body  at  which  the  control  point  is 
located,  and  ACC  is  the  desired  accuracy  of  how  close  a line  source  is  to 
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an  upstream  running  Mach  line  before  it  is  excluded  from  consideration. 
Us\ially,  this  number  is  set  equal  to  approximately  .00001.  Internal  to 
Subroutine  RHP,  the  variable  RHP  is  the  square  of  the  hyperbolic  radius  of 
the  upstream  running  Mach  cone  emanating  from  the  control  point  under  con- 
sideration. 

SUBROUTINE  VELOT  (K,  X,  R,  XO,  RO,  AM,  VX,  VR,  THA) 

Subroutine  VELOT  determines  the  induced  velocity  at  a control  point 
(X,R)  located  on  body  1 due  to  all  applicable  image  line  sources  in  body  2, 
centerline  image  sources  in  body  2,  image  sources  in  body  1 and  centerline 
image  sources  in  body  1.  The  velocity  components  are  summed  and  transferred 
back  as  VX  and  VR  in  the  CALL  argument  list.  The  arrays  XO,  RO,  and  AM  are 
generated  elsewhere  in  the  program  and  are  transferred  into  and  out  of  VELOT 
through  the  CALL  argument  list.  The  variables  X and  R designate  the  axial 
and  radial  coordinate  location  of  the  control  point  under  consideration,  and 
THA  has  the  same  meaning  as  in  Subroutine  RHP.  The  arrays  XO  and  RO  contain 
the  axial  and  radial  location  of  the  line  sources  in  the  upstreams  running 
Mach  cone  from  (X,R,),  and  AM  contains  the  strengths  of  these  line  sources. 

SUBROUTINE  VELl  (K,  X,  R,  SI,  A,  VX,  VR,  THA) 

This  subroutine  determines  the  induced  velocity  at  the  control  point 
(X,R)  on  body  1 due  to  the  distributed  line  sources  (single  body  solution) 
in  body  1.  All  variables  have  been  previously  defined  with  input  information 
entering  through  the  CALL  argument  list . 

SUBROUTINE  VELIT  (K,  X,  R,  SI,  A,  VX,  VR,  THA) 

This  subroutine  is  essentially  the  same  as  VELl  except  it  determines 
the  Induced  velocities  at  the  control  point  (X,R)  on  body  1 due  to  the 
distributed  line  sources  (single  body  solution)  in  body  2, 

SUBROUTINE  VELOCY  (VX,  VR,  TT,  XX,  THA) 

This  subroutine  is  the  sub-executive  program  for  determining  Induced 
velocities  at  the  vatious  control  points . It  calls  the  other  required  sub- 
routines and  transfers  the  summed  velocity  components  and  VX  and  VR  back  to 
the  MAIN  program.  The  control  point  location  under  consideration  is  (XX, TT). 
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REAL  R,L3,LM0,LI,L2.M117,KN 

CI^t£NS^O^)l(C62),R(C60).S(C60)^*ItC60)  ,RT(C6C)»STI060) 

Cl  MENS  ION  CPA(20»,CPB(2C),XE(2C),CP(0t2l 
CCMMO.\/CC/Sn0dO»,A<O6O)  ,RCL2(C60),XIS(060>,RIS1IC60)  , 

*AM2 (060 ) ,XCL (060) |ACL2( Q6C  ) . A K 1 ( C6C ) t RC L I ( 060 ) , ACL  1 (060  ) « 
•RIS2(C60) 

COMMON/AA/N,LMOtCLD 

CONNON/e/Ll ,L2,L3,«MAX, TAN2,A^VIP,M117, VTL 

21  format ( 8 X, IHX, 12X,1HR,13X,1HS) 

22  FORMAT  (1X,3FIA.6) 

29  FORMAT  ( 1 X , 1 3 , AX , AH VR  =F  9.A,AX,AHVX  =F  9.A,AX,AHRC  =F  9. A, 
*AX,3HS  =F9.A, AX, 7HERRCR  =F9 . A , A X , AHCP  = F 1 2 . A , AX  , 3H X =F9.A) 

30  FORMAT ( IX,2H A ( 13, 3H  ) =FIA.6) 

AO  FORMAT ( IX , 15HNC  CF  SOLRCES  =I3,5X,23HNC  CF  PRESSURE  POINTS  = 
♦I3,5X,9HMACH  NC  = F I C . 3 , 5 X , 5HCLC  =F1C.A) 

501  FORMAT  ( 1X,AF2C.5, ICX, 12  ) 

REAC(5, 111)  L1,L2,L3,RMAX,TANG,AAV1P,M117 

111  FORMAT(7F1C.O) 

V«R1TE(6,112)  L 1,L2,L3,RMAX,  TANG,AAV1P,M117 

112  FCRMAT(7F1C.6) 

REAC(5,2)  N,NCP,M,CLD 

2 format  ( 13, I3,F1C.C,F10.C) 

9 REAC(5,7)  THAC, IPLCT,NPLQT 

7 format (FIC.C, 212) 

WRIT£(6,AC)  N,NCP,M,CLD 
Y0=(L1**2-RMAX«*2)/(2.0*RMAX) 

YTL=( (7.17860-L2)**2-RMAX**2)/(2.0*RMAX) 

TA\2=TANG/57, 29578 
LMC=SCRT (N+N-1 .0) 

XLMS=LM0*LM0 
CN=FL0AT  (M 
CX=L3/CN 
X( 1 ) = 0X 
CC  20  J=2,N 
X ( J )=X( j-1 ) iOX 
20  CCNTINLiE 

CALL  6E0ME (R,X,S,N  ) 

V.RI  TE(6,21  ) 

C kRITE(6,22)(X(J),R(J),S(J),J>l,M 

N1=N-1 
SI ( 1 )-0.0 
CC  3 J*2,N 

SI(J)=X( J-1)-LMD*R(J-1) 

C kRITE(6,A9)  J,SI(J) 

C A9  FORMAT  (5X,I2,10X,F1A.6) 

3 CONTINUE 
CO  A I<1,N 
SUN*0.0 
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IF(I.EC.l)  GC  TO  330 
CC  5 

U*(X(  n-SI  (J  ) )/(L^'C♦Rt  I ) ) 

C = ALCG(lJ  + SCRT(L*U-l.O)) 

SL"  =SUR*SI  I )♦#(  J)♦C♦L^'0♦AU)♦SCRT(C♦t-l.C) 

5 CCMINUE 

XNUR  =S(H-SUH 

Ul  = (xn  )-Sl  ( I ) »/(L^'C*R(  M) 

GC  TC  35C 

330  Ul  = X(l)/(LRC*R«m 
XNUR  = R ( I )/X  ( U 
S(  n = XNLK 

350  CENCR  = L^C«SCRT(Ul*tl-l.C)  +S U ) ♦«L0G« L I iSCR  I ( L I *0 1- 1 .0  ) » 
A( I )=XNUR/CENCP 
kR|TE(6.30)  I»A(I) 

A CONTINUE 

C STORE  LOCATICNS  AND  STRENGTHS  CF  IHAGE  SCCRCES  IN  8CCV  I 
C ANC  eCCY  2 AT  ENC  OF  BOGIES. 

CC  50  J=l,N 
XIS(J)=6.C 
R1SI(J)=.0A 
RIS2 ( J)=CLC-.CA 
ARl(J>=C.C 
AR2(J  )=0.0 

C TFE  1 ANC  2 IN  the  variables  ICENTIFIEC  ABOVE  REFER  TO 

C BCCY  1 ANC  2 RESPECTIVELY. 

50  CCNIINUE 

C STORE  LCCATICNS  AND  STRENGTHS  OF  CENTERLINE  SOURCES  IN 

C BCCY  1 ANC  2 AI  THE  ENC  CF  THE  BCCIES. 

CC  51  J=1,N 
XCU  J)=6.0 
RCLl  U)=0.0 
RCL2( J)=CLC 
ACLl(J)=C.C 
ACL2(J)=0.0 

51  CCNIINUE 

IF(CLC.GI.A.O)  GC  TC  l*i 
I = l 

C h£  RUST  NCV.  CETERMINE  WHERE  CN  BCCY  I THE  NCSE  OF  BODY  2 

C FIRST  INTERFERES  WITH  BCCY  1. 

RP=CLC 
XP=C.O 

IF  A PCINT  CN  BOCY  2 IS  GIVEN  IXP.RP)  THIS  SECTION  CF  THE 
PRCGRAR  CALCULATES  THE  INTERSECTICN  CF  THE  CCWNSTREAR 
RUNNING  RACH  LINE  FRCR  IXP.RP)  ANC  THE  UPPER  SURFACE  CF 
BCCY  1.  THIS  INTERSECTICN  PCINT  IS  CENCTED  AS  IXISfRIS) 


i 
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IF  KN  IS  SET  * l.O  , THE  POINT  CN  TCP  CF  THE  BODY  IS 
CCPPLTEC.  IF  KN  IS  SET  » -l.C  t THE  PCINT  CN  THE  BCTTCP 
CF  THE  eCCY  IS  CCPPLTEC. 

KN=  l.C 

XINaXP*(Ll/RP#X»*(  {LMC*CLCU(LPC+(Ll/ftPAXn  » 
kRITE(6.1l3)  XIN 
113  FCRPAT(AHXIN=,F10.6) 

70  KCNT=1 
Y*X1N 

eo  IFIXIN.lt. Ll)  RIN=-YC»SCRT( (L1**2+YC**2)-«L1-XC I ) )**2) 
IF( X IN.GT .Ll .ANO.X IN.LT .L2)  KIN^RPAX 
IF(XIN.GT.L2)  GC  TC  1A 
XN£V.  = XP*lPC*(RP-KN*RIN) 

IF(AeS(XlN-XN£h).LT.CCCCCl)  GC  TC  210 

C IF2  = X IN-XNEH 

IF(KCNT.£C.2)  GO  TC  190 

CIF1=CIF2 

XIN=XN£W 

KENT  = KCNT+l 

GC  TC  0C 

19C  XIN=(CIF1*XIN-CIF2*Y)/(0IF1-0IF2I 
GC  TC  70 

210  IF(KN.LT.C.CIRIN=-RIN 
XIN^XNEki 
NAP=1 

V«RITE(6tS01  ) XP,RP,X1N,RIN.NAP 
CC  60  J=ltN 

IFI  XIN.GT.XUn  GC  TO  60 

JSAVE  » J 

XTSAVE=X(J» 

RTSAVE-RI J» 

SSAVE»SU) 

WRITE (6. 69)  JSAVE, XT  SAVE. RTSAVE.SSAVE 
69  FCBPAH5X,I2,5X,3FIA.6) 

GC  TC  62 
60  CONTINUE 
62  XT(n*XTSAVE 
RTI  I )<RTSAVE 
ST(  DISSAVE 

RISl  1 1 )».5*MXTSAVE-XS)/LPD*(RS-RTSAVE)  ) 

KISI  I )*XS-LPC*(RS-RISU  I) ) 

XCU  I )*X  H I l-LPOARTI  I ) 

RIS2m^CLL-RISim 

IP«JSAVE«1 

kRITE(6,S0S) 

SOS  FCRPAT(19X,3HXIS, 17X,3HRIS,17X,3HXCL,12X,1HI) 

MRITEI6,SCA  IXISM  ),RIS1  ( n.XCLI  I 1,1 


5C4  FCR^AT( IX.3F20.6, 110) 

5 I2>l«l 

XTU2)  = X»  IP) 

RTI  I2)>RI IP) 

STII2)*S( IP) 

RISUI2)=.5*((XT(I2)-XTm)/LPC*IBT(I)-RTn2))) 
XIS(I2)=XT(I )-LMC*(RTI I )-RISll 12) ) 

XCL(I2>=XT( 12 )-LPC*RTI12) 

RIS2I  12)=CLC-RIS1M2) 

WRITE(6,5C4)XIS(l2).RISl(12lfXCL(I2).I2 
IF( xri 12) .GT.5.99)  GO  TG  56 
1 = 1 + 1 
IP=IP+1 
GC  TC  55 
56  JAP=I 

kRITE(6t66  ) 

6 6 FCRP4T  ( 12X,  2HXf',  12X  ,2HRP,  12X,  3hXhH,  11X,2HXX,  12X,2HT  r,12X, 
*3F6Pl,llX,2FXe,l2X,2HeBf IlX.AFACL  U 
CC  ICC  J=  1,  JAP 
VX  = C.O 
VR  = 0.0 
XX=XT( J) 

TT=RTIJ) 

XB=XT(J+1) 

EE=  RT  ( J + 1 ) 

IFITT.LT.RISl  (jn  GC  TC  ICO 
CALL  VELCCV  (VX,VR,TT,XX,1.57C796) 

IF(AeS(VX).LT..CCCCCl)  GC  TO  ICl 

NCTE  AT  This  PCINT  J IS  TEE  NOPfifcR  CF  TEE  CCNTRCL  PCINT 
T)-AT  IS  FIRST  AFFECTED  BY  THE  PRESENCE  CF  THE  OTHER  8C0Y. 
XP  = XX-XISU) 

HP=RIS2IJ)-TT 
XHHsXP+XP-LPC+LPC+RP^RM 
XNUP=ST(J)*I l.C+VX)-VR 

ag=xp/(lpc*itt-risi(j))) 

SRT=SCRT(AG+AG-1.0I 
IF!  XHH63,63,6^ 

6 3 APHJ  ) = XNLP/ILPD*SR1+STI  J)*ALCG(AG+SRT)  ) 

GC  TC  65 

6<*  AGG  = xr'/(LPC*(RrS2lJ)-TTI) 

SRTT=SCRI(ACG*AGG-1.0) 

RCX=LPC*ISRT-SRTT ) 

RCC=ALCGI (ACG+SRTT)*IAG+SRT)) 

API! J)=XNLP/I RCN+ST (J)*RCC) 

65  AP21J)=AP1 (J) 

RE  NCR  CCNSrCER  THE  CCNTRCL  PCINT  ON  THE  OCTTCM  OF  BODY  I 
01  VX=C.C 
VR=C.C 


X’ 


s 

I 

i 

[ 


[ 

i 


i 


CALL  VELCCV (VXtVHtee.XB  fA. 712369) 
|F(A6SIVX).LT..CCC001  ) GC  TO  ICC 
AGB*(xe-xcujn/aHC*ee) 

SRIBaSC(»T(4GP*AGH-l.O) 

XNUP»ST) J*1 )*) l.0*VX)-VR 

CENO'  = LPC*SRTe*Sn  J*l  )*ALCC(  AGH+SRTB  ) 

ACL1U)»X^LP/CE^C^» 

ACL2IJ)=ACL1IJ) 

hRITE  (6«67  ) XP,RH,Xt-H,XX,TT.ARl(  J),XB,BP,ACLl(J) 
67  FCRPAT«2X,9F1«.5) 

ICC  CCMINUE 
14  C\N»FLCAT<NCP) 

CX=L3/CNN 
X( 1 |«CX 
CC  6 J*2iNCP 
XJJ)*X(J-1)*CX 
6 CCMINLE 

CALL  GECRE«R,X,S,NCP) 

WRITE  (6.21) 

WRl  IE(6.22  ) (X(J).R(J)  tSU).J>l,^CP) 

IF(  IPLCT.EC.O)  GC  TC  105 


CALL  GSI7E(e.C.ll.C.llC0) 
CALL  PLCn2.5,7.0,-3) 


1C5  KS=C 

TFAC=THAC-45.0 
CC  16  J»l,NPLOT 
TFAC  = TFAC+45.0 
WRITE(6,P)  THAC 
e FCKPAT(20X,7HfFETA  *FI0.2> 

TFA=TFAC/57. 29570 

CC  13  1*1, ^CP 

VX=C.C 

VR=C.O 

TT  = R(  I ) 

XX=X(I) 

CALL  VELCCY  I VX , VR . T T , XX , THA ) 
eC=VR/( l.C»VX) 

ERRCR=(PC-SI I ) )*57. 29578 
CPI n*-2.0«VX«LPC*LPC*VX*VX  -VRAVR 
IFICPI  I ) .GT.  2.0)  CPI  I )«  2.0 
IFICPI!  ).LT.-2.0)  CPII)«-2.0 
WRITE  16.29)  I.VR.VX.BC.SI D.ERRCR.CPI l).X(I) 
13  CCMINLE 

IFI  IPLCT.EC.O)  CC  TC  16 
CC  17!=l,^CP 
CPI  I )*CPI  I lAlC.O 
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2 <1 

\ 1 

1 

1 


17  CCMINUE 

CC  le  I»I,15 
CPA( I )=CP«( 1 )*10.0 

cpet  n=cpe(  n*io.o 

18  CCMINUE 
NCPl=NCP*l 
NCP2=NCP*2 
CPINCPl  )-C.O 
CP(NCP2)=1.0 
X(NCP1)=0.0 
X(NCP2)=1.C 
RINCPl)=O.C 
R(NCP2)=l.C 
CALL  FACTCRI.eO) 

CALL  AxlSIC.C,C.0,lHX,-l,6.0,C.C,XlNCPl»,XINCP2n 
CALL  AX IS(C.C,-2.0, 5HCPX1C,5,5.C,9C.C,-2.C,CPINCP2»  I 
CALL  SYPeCL (2.C.2.5.  . lC,2ChA9V  IP-M= 1 . 5 1- ALP  = C .0 ,0 . C t20  ) 
CALL  SYPeCL(3.Ct2.Cf.  1C,  lAHEXPERIPENrAL"tC.CfI«) 

CALL  SYP0CL(99S.O,2.O5,.1C,C,C.C,O) 

CALL  SYMDCL  (3.C,  1.7, . lC,tHHECR>,0.0,6) 
call  LINEIX, CP, NCP,  1,0.0) 

CALL  PL0T(C.C,-9.C,-3) 

CALL  axlSIC.CfC.Ct  U*X,-1  ,6.CrC.C,XlNCPl) ,X(NCP2)) 

CALL  AXISIC.C,C.0,lHR,l,l.C,9C.C,R(NCPl),RINCP2n 

x(  n=c.o 

R(  1 ) = C.O 

CALL  LINE(X,R,N, 1,0,0) 

CALL  PLCT(E.0,C.0,2) 

CC  129  IK=1,N 
RI  IK)=-RI  IK) 

129  CCMINUE 

CALL  PLCnC.C,C.0,3) 

CALL  LINE(X,R,N,  1,0,0) 

CALL  PLCT(6.0,C.0,2) 

CALL  PLCKC.C.C. 0,999) 

16  CCNTINUE 
STCP 
ENC 


r ' 


SUBRCLTINE  GECf»e  ( RAC,  X,  SLOPE  tNS  ) 

REAL  LltL2,L1tPll7 
CIP ENSIGN  RACiNSItXINS) .SLOPE  I NS  I 
CCf'PCN/e/Ll,L2,L3.RPAX,TAN2,AAWlP,Kll7,VTL 
CC  SSq  I*  l.NS 
IF(X(  n.GT.Ll)  GC  TO  ICCC 

: NOSE  SECriCN 

YC  = IL1**2-RPAX«*2)/ I2.0ARKAXI 
RACm=-VC*SCRTMLl**2»YC**2)-(Ll-X(  I n**2) 

SLCPEII  ) = (Ll-X(  I ) )/SQRTMLlA*2*YC*A2)-(Ll-Xm)**2) 

GC  TO  qqq 
: PIC-SECTICN 

ICCC  IFIXI  I).GT.L2»  CC  TC  2CCC 

RAC  ( I I^KPAX 
SLOPE!  n=c.c 
CC  TC  qqq 
: TAIL  SECTION 

2CCO  IF!X(  1I.GT.L3I  CC  TC  888 

RAOn)  = (-YTL  + SCRT((  ( 7 . 1 7 860-L  2 I * *2  ♦ Y TL*  *2  ) - ( X 1 1 ) -L2  I ♦♦2 ))  • 
♦AqvIP*  (KPAX-IAM2*! X(I)-L2))*PII7 
SLOPE (I ) = (-(X( n-L2  )/SCRT( ( ( 7 . I TE60-L2 I ♦ YTL**2 I - ( X (I )- 
♦L2)**2) )*AAVIP  - (TAN2*P117) 

GC  TC  qqq 
EEP  RAC  ( I ) = C .C 
SLCPEI I )=C.O 
ssq  CONTINUE 
RETURN 
ENC 


o n o r» 


St‘fiRCLTl^E  RhP(K,X,Q,XI,RI,THA,«CC) 

REAL  LPCtKN 
CCRRCN/AA/K,LPC,CLD 
CIPENSICN  XIIN),RIIN) 

Al  FCRPATUX.IEK  *I2» 

TUS  Sl;eRCtTt^E  CETERMINES  HOh  AA\V  OF  TFE  SfACULARIMES 
IK  8CCV  1 CR  2 IKFLtENCE  ThE  CCMRCL  PCINT  XIJI,  R(J)  CK 
eCCV  1.  TFIS  KUPRER  IS  STOREC  CR  TRANSFERRED  AS  K. 

V * -R*CCS(THA) 

2 * RPSINITFA) 

CC  lA  JJ^UK 
IFIXIUJI.GT.XIGC  rc  18 

HRP>|X-X|  I jjn*«2-LFOPLPCAI  (O.C-V)*«2«IRIIJJI-n*«2)-ACC 
IFIFRP.lt. C.O)  GC  TC  18 
K = K*l 

1<>  CCKTIKUE 
C IS  FRITElEtKlI  K 
18  HETLRK 
EKC 
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SUPBCCTINE  VELCnK,X,B,XC,KC,*H,VX,VR,ThA) 
REAL  IPC 

CCRRCN/AA/N,LPC,CLC 
tIPENSICN  XC(N)«RO(N).AHIN) 

IF(K.EC.O)  GC  TO  18 
CC  16  J^ltK 

Cl)»*V  * RCIJ  )*RCIJ»*R*B 

RP  - SCRHCLPY-2.0*B*PQI  J)*SIN(THAn 

AG  *l*-XCIJ))/(LPC*BP) 

SRT=SC«r(AG*AG-1.01 

VX»VX-APtj)*ALCG(AG*SRr) 

VRP*  AP(J  l*LPC*SRT 

si\e=R*ccsnFA)/Rp 

CCSB*(RCIJ)-R*SIN(THA))/RP 

vz=-vpp*ccse 

VY=-VPP*S lAB 

VR=VR*VZ*S iNCThA )-VY*COS< THAI 
16  CCNTINGE 
18  RETLRA 


SUBRCLT INE  VEl U K , X ,R, S I , A, VX , VR t THA ) 
real 

CC^'^'C^/AA/^,L(»C,CLC 
tlHEKSlCN  SnN),A(N) 

IF!K.EC.0»  GC  to  17 
CC  16  J=1,K 
AG=  (X-St(J))/(Lf'C*R» 
SRT=SCRT(AG*AG-1.0) 

VX*VX-A( J)*ALCGrAG*SRTI 
VR  = VR+A  < J )'*L^'C*SRT 

16  CCMINLE 

17  RETURN 
ENC 


J 

j . SLBRCLTINE  VEL I T ( K, X ,R, S I , A , VX , VR, THA » 

REAL  LMC 

CC^<^'C^/AA/^,L^'C,CLC 
CIRE^SIC^  SI(N),A(N) 

’ IFCK.EC.C)  GC  TO  17 

VRP=0.0 

CUPY=CLC*CLC+R*R 

RP=SCRT ICLPY-2.0*R*CLC*SIN( IHA) I 

CC  16  J=l,K  I 

AG=  U-Snj))/(LPC*flP>  ! 

SRT=SCRT <AG*AG-l.O»  i 

! VX=VX-A( J)«4LCG<AG*SRT)  ■ 

j VRP  = VRP  + A (J  )♦L^'C♦SRT 

' 16  CCNTINUE  i 

SINP=R*CCS(TBA)/RP  ] 

CCSe=(CLC-R♦£I^(T^•A))/RP  ! 

I V2=-VRP*CCSe  ? 

VY=-VRP*SINe  j 

VR=VR+VZ*S IN( ThA )-VY*CCS( THAI 

17  RETLRN  j 

ENC  i 


I 


I 


i 


I 


1 

I 


SCRPCLTINE  VELCCY  ( VX ,VR , T T , XX , THA » 
CCKr>CN/CC/SMC6C)tA(0£0)«RCL2(C(C).XIS(C6Ci.RlSllC6CI  t 
• A><2(C6C).XCUC60).ACL2(C6C)f  AM(C60>|RCLI(C60),ACL1(C6CI. 
•RIS2IC60) 

C KLST  flRSI  CETtRf'CNE  HGh  RASY  LINE  SOLRCE  SEGMENTS  INFLLENCE 

C CHECK  INCUCEC  VELOCITIES  FRCR  LINE  SCLRCES  IN  BCCY  2.  ViE 

C THE  CCNTrtCL  PCINT  J 

K1  > 0 

call  RHP  (Kl,  XX,  TT,  SI,  RCL 2 , THA , .CCO CC 1) 

CALL  VELITIKl,  XX,  TT,  SI,  A,  VX,  VR,THA) 

C ViE  HLST  NCV.  CETERPINE  HOW  PANV  LINE  SCLRCES  ALCNG  THE 

C IPACE  LINE  CF  BCCY  2 INFLLENCE  THE  CCNTRCL  PCINT  XU),  R(J) 

C NCW  LNCER  CCSSICERATICN. 

K2  * C 

CALL  RHP  (K2,XX,TT,XIS,RIS2,THA,.CCCCCl ) 

C WE  PLST  NCW  CETERPINE  THE  INCLCEC  VELCCITIES  FRCP  THE  LINE 
C SCLRCE  ALCNG  IHE  IPAGE  LINE  IN  ECOY  2. 

CALL  VELCT IK2,XX,TT  ,XIS,RIS2, AP2,VX,VR, THA) 

C WE  PLST  NCW  CETERPINE  HCW  PANY  LINE  SCLRCES  ALCNG  IHE 

C CENTERLINE  CF  BOCY  2 INFLLENCE  IHE  CCNTRCL  PCINT  XU),  RIJ). 

K3  = 0 

CALL  RHP  (K3,XX,TT,XCL,RCL2,THA,.CCCCCI) 

C WE  PLST  NCW  CETERPINE  THE  INCLCEC  VELCCITIES  FRCP  LINE 

C SCLRCES  ALCNG  IHE  CENTERLINE  CF  BCDY  2. 

CALL  VELCT(K3,XX,TT,XCL,RCL2,  AC L 2 , VX , V R , TH A ) 

C WE  PLST  NCW  CETERPINE  THE  INCLCEC  VELCCITIES  AT  THE  CCNTRCL 

C PCINT  CLE  THE  CISTRIflLTEC  SCLRCES  IN  BCCY  I . 

K6  * 0 

call  RHP(K6,XX,TT,S  I,RCLl,THA,.CCCCOl) 
call  VELllKe,XX,TT,SI,A,VX,VR,THAI 
C WE  PLST  NCW  CETERPINE  THE  INCLCEC  VELCCITIES  AT  THE  SAPE 
C CCNTRCL  PCINT  FRCP  THE  IPAGE  LINE  ANC  CENTERLINE  LINE 
C SCLRCES  IN  BOCY  1. 

KA  = C 

CALL  RHP  IKA,XX,TT,XIS,RISI,THA,.CCCCC1 ) 

CALL  VELCT (KA,XX,TI ,XIS,RISI,API,VX,VR, THA) 

K5  =0 

CALL  RHP  (K5,XX,TT,XCL,RCLl,THA,.CCCCCl) 

CALL  VELCT (KS, XX, TT  ,XCL, RCL  1,  ACL  1 , V X , VR , TH A ) 

RETLRN 

ENC 
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